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Abstract. Wavefront sampling experiments record arrival 
times of Cerenkov photons with high precision at various lo- 
cations in Cerenkov pool using a distributed array of tele- 
scopes. It was shown earlier that this photon front can be 
fitted with a spherical surface traveling at a speed of light 
and originating from a single point on the shower axis. Ra- 
dius of curvature of the spherical shower front (R) is approx- 
imately equal to the height of shower maximum from obser- 
vation level. For a given primary species, it is also found 
that R varies with the primary energy {E) and this provides a 
method of estimating the primary energy. In general, one can 
estimate the arrival times at each telescope using the radius 
of curvature, arrival direction of the primary and the core 
location. This, when compared with the data enables us to 
estimate the above parameters for each shower. This method 
of obtaining the arrival direction alleviates the difficulty in 
the form of systematics arising out of the plane wavefront 
approximation for the Cerenkov front. Another outstanding 
problem in the field of atmospheric Cerenkov technique is the 
difficulty in locating the shower core. This method seems to 
solve both these problems and provides an elegant method 
to determine the arrival direction as well as the core loca- 
tion from timing information alone. In addition, using the 
Cerenkov photon density information and the core position 
we can estimate the energy of the primary if the nature of 
the primary is known. Combining these two independent es- 
timates of the primary energy, the energy resolution can be 
further improved. Application of this methodology to sim- 
ulated data and the results will be presented. The intrinsic 
uncertainties on the various estimated parameters also will 
be discussed. 



experiment. The density measurements and core location in- 
formation enable one to estimate the primary energy which in 
turn could be used to derive the energy spectrum of primary. 
On the other hand, the timing information which enables us 
to derive the arrival direction of the primary is essential for 
source search. 

Estimation of shower parameters is rather easily done at 
PeV energies by measuring particle densities and their ar- 
rival times. However at TeV energies, where the charged 
particles do not reach the observation level, the measurement 
of shower parameters is rather a challenging task. This is 
primarily because a bulk of the recorded events have the core 
locations far away from the points of measurement. In the 
case of 7-ray primaries, the Cerenkov photon lateral distri- 
bution does not exhibit significant gradient which makes it 
even more difficult to estimate the core location. Further, in 
the case of Cerenkov telescope arrays, in most cases, only a 
part of the Cerenkov pool is intercepted by the array. As a 
result the estimate of the energy of the primary will be very 
uncertain without the core location information. 

The arrival direction of the primary can be estimated from 
the timing measurements of the Cerenkov light front at spaced 
telescopes. The conventional plane front approximation leads 
to significant systematic errors in the estimated angles. These 
errors often undermine the otherwise excellent angular reso- 
lution of the system. In order to eliminate the interference 
due to systematic errors one has to use a spherical approxi- 
mation to the light front. As a result, we can obtain impor- 
tant additional information like the radius of curvature of the 
Cerenkov front and the core location. 

Here we present the methodology and results of such an 
effort based on simulation studies. 



1 Introduction 

Measurement of shower parameters like the primary energy, 
core location and direction of arrival are vital in any air shower 
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2 Cerenkov shower front fitting for vertical showers 

In a typical wavefront sampUng experiment arrival time of 
Cerenkov shower front is recorded at several locations in 
Cerenkov pool with high precision. This information can be 
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Fig. 1. Variation of mean arrival time of Cerenkov shower front 
as a function of core distance for 7— rays, protons and Fe nuclei of 
various energies incident vertically at the top of the atmosphere. For 
each case, arrival time is averaged over typically 10 showers. The 
smooth curve corresponds to the best fit spherical wavefront. 



used to reconstruct the shower direction. It is shown earlier 
that for vertically incident showers initiated by 7— rays and 
cosmic ray primaries, the relative arrival time delay [t{r)] of 
Cerenkov shower front at a core distance r can be approxi- 
mated by 



V(i?2 +r2) R 



(1) 



where R is the radius of curvature of the spherical front (Bat- 
tistoni et al. (1998), Chitnis and Bhat (1999)). 

In effect, Cerenkov shower front can be approximated with 
a wavefront moving at a speed of light, originating from a 
single point on shower axis. Figure 1 shows the variation 
of mean arrival time of Cerenkov shower front as a function 
of core distance for 7— rays, protons and Fe nuclei of vari- 
ous energies incident vertically at the top of the atmosphere. 
These showers were simulated using a package known as 
CORSIKA (Heck et al, 1998). For each case, arrival time is 
averaged over typically 10 showers. Best fit spherical wave- 
front is also indicated by smooth curve. Radius of the fit- 
ted wavefront (R) approximately corresponds to the height 
of shower maximum (he) from observation level. This is ex- 
pected since large fraction of emission arises from the region 
in the vicinity of shower maximum. 

3 Algorithm for spherical wavefront fitting for inclined 
showers 

Here we extend this approach to the showers generated by 
primaries incident at an angle with respect to the vertical. 
Consider a distributed array of telescopes at observation level. 
Shower is incident at a zenith angle 9 and azimuth angle <j). 
Here it is assumed that the shower axis passes through the 
centre of the array. Coordinates of the telescopes in the array 
are given by ixi,yi), with respect to the centre of the array. 



The arrival time of the shower front at a telescope situated 
at (xi,yi), with respect to that at the footprint of the front, is 
given by 



tir) = 1 



i?2 



- 2Rtan6 [xicoscf) + yisincp) +yf — R 



(2) 



where R is the vertical height of the emission point on shower 
axis from the observation level. 

Using this algorithm one can fit the measured arrival time 
and estimate the arrival direction of shower front (6 and (p) 
as well as the radius of curvature of the shower front (R'). 

4 Spherical wavefront fitting of simulated data 

Efficacy of this algorithm is tested using simulated data. A 
large number of showers initiated by protons incident at an 
average angle of 15° with respect to vertical are simulated 
using CORSIKA (Heck et al., 1998). Energies of these pro- 
ton primaries are chosen randomly in the range 250 GeV - 
20 TeV with the spectral index of -2.65. Zenith angles for 
showers are chosen randomly within 2° around the average 
value. Whereas azimuthal angles are randomly selected in 
the range 0-360°. Axes of all the showers pass through the 
central telescope of the array. Arrival times of Cerenkov pho- 
tons are recorded at 357 telescopes spread over an area of 400 
m X 400 m at an altitude of 1 km. 

For each shower, arrival times at various core distances are 
fitted with a spherical wavefront. Angles 9 and (j) and verti- 
cal height of emission point from observation level (i?) are 
estimated using Marquardt routine. This exercise is carried 
out for a sample of 50 showers. Figure 2(a) shows the dis- 
tribution of fitted radii of curvature (R) of shower front. For 
comparison, distribution of the shower maxima hg from ob- 
servation level is also shown in the same figure. Distribution 
of difference in R and hp is shown in Figure 2(b). Mean 
value of this difference distribution is 3.3 km with an RMS 
spread of 2.8 km, indicating that effective point of emission 
is above the shower maximum. Difference distribution of 
fitted azimuthal angle ((j)) and the actual azimuthal angle is 
shown in Figure 3a. Similarly difference distribution of fitted 
zenith angle (9) and actual zenith angle is shown in Figure 3b. 

So far we have carried out spherical wavefront fitting as- 
suming shower core to be exactly at the centre of the array, 
which is also the origin of our coordinate system. However, 
in an actual experiment shower core may lie anywhere in the 
telescope array and at times even outside the array. By intro- 
ducing the core location as another variable, same algorithm 
can be further extended to get location of the shower core. 

5 Discussion and conclusions 

Cerenkov shower fronts from proton showers incident at an 
angle of 15° are fitted with a spherical wavefront. Fitted ra- 
dius R is found to be less than the height of shower maximum 
he by about 3 kms. For vertical showers R is found to be 
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Fig. 2. (a) Distribution of tlie fitted radii of shower front R (con- 
tinuous fine) and distribution of the actual shower maxima he. from 
observation level (shown by dotted line), (b) Difference distribution 
of R and fte. Mean of the distribution and the expected value are 
marked by dashed and dotted lines. Mean and RMS of the differ- 
ence distribution are indicated. 



about 1 km lower than he (Chitnis and Bhat, 2001). Larger 
difference between R and he for inclined showers could be 
intrinsic. Reduction in the relative contribution of the cur- 
vature to the arrival time differences at various telescopes 
compared to those due to propagation delay could result in a 
larger error in fitted radius as well. One can subtract the rel- 
ative time of arrival differences arising purely out of spatial 
separation of the detectors using plane front approximation. 
The residuals can then be fitted to a spherical front which 
could enhance the sensitivity of the curvature of the shower 
front which in turn will reduce the systematic error on the 
fitted radius of curvature. We plan to implement this in our 
future fitting procedures when we will be able to present a 
quantitative estimate of the degree of improvement in the fit- 
ted radius of curvature. 

For a given species the height of shower maximum de- 
creases as the logarithm of the primary energy (Rahman et 
ai, 2001). This offers us an alternate method of estimating 
the primary energy which in turn improves the primary en- 
ergy estimate. 

In addition to the above, spherical fit to the Cerenkov light 
front offers arrival direction as well as core location infor- 
mation of the shower As mentioned in the introduction the 
arrival direction information derived by this method is rela- 
tively free of systematic errors and hence one can fully ex- 
ploit the actual angular resolution of the experiment limited 
by the system hardware. This, in turn offers a major advan- 
tage of rejecting a bulk of the off-axis showers which could 
be presumed to be due to cosmic rays with better confidence. 

Finally, shower parameters estimated by this methodol- 
ogy like the radius of curvature and core location coupled 
with measured and simulated photon density distributions 
will greatly improve the estimation of the primary energy. 
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Fig. 3. (a) Difference distribution of fitted and actual Mean 
of the distribution and the expected value are marked by dashed 
and dotted lines. Mean and RMS of the difference distribution are 
indicated, (b) Difference distribution of fitted and actual 9. Mean 
of the distribution and the expected value are marked by dashed 
and dotted lines. Mean and RMS of the difference distribution are 
indicated. 
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